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ABSTRACT 

We examined the variability of three ultra- himinous X-ray sources (ULXs) in the 2003, 
110 ks XMM-Newton observation of NGC 253. Remarkably, we discovered ULXl to 
be three times more variable than ULX2 in the 0.3-10 keV band, even though ULX2 
is brighter. Indeed, ULXl exhibits a power density spectrum that is consistent with 
the canonical high state or very high/steep power law states, but not the canonical 
low state. The 0.3-10 keV emission of ULXl is predominantly non-thermal, and may 
be related to the very high state. We also fitted the ULX spectra with disc blackbody, 
slim disc and convolution Comptonization (siMPL(g)DiSKBB) models. The brightest 
ULX spectra are usually described by a two emission components (disc blackbody 
-f Comptonized component); however, the SIIVIFL model results in a single emission 
component, and may help determine whether the well known soft excess is a feature 
of ULX spectra or an artifact of the two-component model. The SIMPL models were 
rejected for ULXS (and also for the black hole -1- Wolf-Rayet binary ICIO X-1); hence, 
we infer that the observed soft-excesses are genuine features of ULX emission spectra. 
We use an extended corona scenario to explain the soft excess seen in all the highest 
quality ULX spectra, and provide a mechanism for stellar mass black holes to exhibit 
super-Eddington luminosities while remaining locally sub-Eddington. 
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1 INTRODUCTION 



In lBarnard et al.l (|2008al ') we identified a new spectral state 
in the confirmed black hole -I- Wolf-Rayet binary ICIO X-1, 
as well as the BH-I-WR candidate NGC 300 X-1 and the 
weU known black hole binary LMC X-1. Both ICIO X-1 
and NGC300 X-1 exhibited ~ 0.0001-0.1 Hz power den- 
sity spectrum (PDS) where the power is inversely propor- 
tional to the frequency, like the canonical high state. How- 
ever, the 0.3-10 keV for these sources emission is predomi- 
nantly non-thermal, while the canonical emission spectrum 
for high state black hole X-ray binaries is dominated by a 
multi-temperature disk blackbody that contribute s ~90% of 
the 1-10 keV flux ()McClintock fc RemiUardll2006h . 

In iBarnard et all l|2008ah we proposed that the non- 
thermal high state was due to a persistently high accretion 
rate maintaining a stable corona; most black hole X-ray bi- 
naries are transient, and we suggested that the canonical 
soft high state was caused by the eje ction of the corona dur - 
ing outburst. After finding out from ISoria fc Kuncid (|2008l ) 
that ~90% of X-ray sources with luminosities >10''^ erg s~^ 
have non-thermal spectra, we suggested that the same phys- 



ical processes may apply in these sources as in ICIO X-1 and 
NGC 300 X-1. The ultra-luminous X-ray sources (ULXs) are 
an intriguing class of extragalactic X-ray sources that are not 
associated with the galaxy nucleus, yet exhibit X-ray lumi- 
nosities >2xl0''^ erg s — 1, the Eddin gton limit for a ~15 
Mq black hole (see e.g. lRobertsll2007l . for a review). 



Our survey of X-ray sources in the 2003 XMM- 
Newton observation o f the nearby spiral galaxy NGC 253 
(jBarnard et all l2008bh yielded three X-ray sources with 
0.3-10 keV luminosities >2xl0^^ erg s"\ We shall caU 
these ULX1-ULX3, and their known properties are pre- 
sented in Table l T] thes e properties are taken from our survey 
(|Barnard et allboOSbl ). In this paper, we examine the vari- 
ability of these sources. We also perform additional spectral 
fitting; we model their spectra with slim disc models, and 
"physical" Comptonization models where the inner disc pro- 
vides the seed photons. In Section 2 we discuss the observa- 
tion and data analysis, followed by the results in Section 3. 
Our discussion and conclusion is presented in Section 4. 
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1.1 ULX variability 

Few ultra-luminous X-ray sources (ULXs) have sufficiently 
high quality data for timing analysis. The PDS of M82 X- 
1 has been characte rised by a quas i-periodic oscillations 
(QPOs) at 54 mHz jStrohmaver fc M ushotzkv 200^) and 
114 Hz (|Dewangan et a l. 2006b), on top of a broken power 



law continuum ( Dewangan et al]|2006bl ). QPOs have also 



been fo und in Holmberg IX X-1 (along with a power law con- 
tinuum [DewmganeEaO[2006i}_andjre a 



bro ken power l a w con tinuum Strohmaver et al]l2007l ) 



iHeil et all (|2009l ) conducted a timing study of ULXs 
observed with XMM-Newton. They analysed 19 observa- 
tions of 16 ULXs and found 6 ULXs with PDS best de- 
scribed by power laws or broken power laws: M82 X-1, NGC 
1313X-1, NGC 1313 X-2, NGC 55 ULX, Ho IX X-1 and 
NGC 5408 X-1; t he others showed no significant variability. 
iHeil et al.l l|2009l ) obtained upper limits to the variability of 
these ULXs, and concluded that the intrinsic variability was 
lower. 



1.2 ULX emission spectra 



Two emission components are often used to describe the X- 
ray spectra of ULXs: a thermal component, and a compo- 
nent, e.g. power law, to represent inverse C ompton scatter- 
ing o f cool photons on hot ele ctrons, (see e.g. lStobbart et eiD 
I2OO6I : [Gladstone et al.l 120091 '). These fits often result in a 
"soft excess" , where the power law component diverges from 
the thermal component Such spectra have been deemed 
by some as non-standard and unphysical, assuming that 
the seed photons are provided by the inner disc (see e.g. 
iRoberts et al.ll2005l: IConcalves fc Soriall2006t ). 

Gladstone et al I (|2009l ') conducted a survey of XMM- 



Newton observations of ULXs with > 10000 source counts 
in their EPIC spectra; their preferred model consists of a 
disc blackbody and a Comptonized component with the seed 
photons tied to the inner disc temperature. They found two 
universal features: a soft excess and roll-over in the spectrum 
above ~3 keV. They found that disc + Comptonized corona 
emission models fitted the data well; however, the resulting 
coronae were cool (kT ~ 1-3 keV) and optically thick (r ~ 
5-30). The combination of a cool d isc and broken hard com- 
ponent led [Gladstone et alj |2009l ) to define a new, ultralu- 
minous spectral state. Such coronae are distinct from coro - 
nae of Galactic binaries (r ~ 1), and lGladstone et al.l (|2009l ) 
use this difference to reject the scenario of sub-Eddington 
accretion onto intermediate mass black holes. 



2 OBSERVATION AND ANALYSIS 



In lBarnard et all (|2008bl ) we conducted a spectral survey of 
X-ray sources in the 2003 XMM-Newton observation of NGC 
253; full details of our data reduction are given in that paper. 
The prese nt work concerns new variability studies of three 
ULXs. In 'Barnar d et all (|200it ) we showed that the combi- 
nation of lightcurves from the different cameras on-board 
XMM-Newton can lead to artefacts in the PDS. To ensure 
fidelity of the PDS, we performed our variability analysis 
on the unbinned pn events files for each source. The total 
exposure time was ~110 ks; however, several background 



Figure 1. Detail of a three colour pn -|- MOS image of NGC 253 
showing ULXl, ULX2 and ULXS. Red represents 0.3-2.0 keV, 
green represents 2.0-3.0 keV, and blue represents 4.0-10 keV. 
The image is log-scaled, and the source extraction regions for the 
ULXs are labeled. 



flares occur red during the observ ation; the flares were re- 
moved (see iBarnard et al.|[2008bl . for details), leaving ~50 
ks of good data. 



3 RESULTS 

Figure [1] shows a three colour detail of the combined EPIC 
(pn + MOSl -f MOS2) image of NGC 253, overlaid with the 
source extraction regions; the image is log-scaled, north is 
up and east is left. 



3.1 Variability 

We present the 0.3-10 keV pn lightcurves of ULXl, ULX2 
and ULX3 in Figure (2] binned to 100 s from the unbinned 
events. Intervals with strong background flares have been 
flltered out. The background contribution is negligible in 
each case. The x- and y-axes are plotted on the same scale 
for each source. For each lightcurve we provide the fractional 
r.m.s. variability calculated by the program lcurve from 
the FTOOLs analysis suit^B- We note that ULXS is rather 
fainter than ULXl and ULX2 and there are two reasons 
for this. Firstly, ULXS is located at the corner of a chip in 
the pn image, and the collecting area is reduced when the 



http: / /heasarc. nasa.gov/lheasoft/ftools/ftools_menu. html 
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Table 1. Known properties of the three ULXs, taken from our survey of X-ray sources in NGC253 l|Barnard et al.ll2008d ). For each 
source we give the source coordinates as returned by the source detection routine; these have uncertainties of ~2". We then give the 
number of net source counts in pn and MOS detectors. Next we give the details of the best fit spectral model: absorption (A'^h), blackbody 
temperature (kT) and photon index (F), along with the corresponding x^/dof and 0.3—10 keV luminosity / 10^^ erg s~^. Numbers in 
parentheses indicate 90% confidence limits on the final digit. 



Source 


Position 


pn 


MOS 


JVh/10^^ atom cm"^ 


kT/keV 


F 


xVdof 


L/1039 


ULXl 
ULX2 
ULX3 


00 47 22.56 -25 20 51 
00 47 32.98 -25 17 50 
00 47 42.41 -25 14 59 


10537 
11614 
773 


10390 
12014 
4156 


2.00(2) 
2.9(2) 
6.0(11) 


0.73(5) 
0.98(6) 
0.94(8) 


2.14(4) 
1.94(5) 
3.4(5) 


347/323 
374/374 
84/76 


2.90(12) 
4.10(19) 
2.4(4) 
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Figure 2. Lightcurves of the three ULXs, binned to 100 s from 
the unbinned 0.3—10 keV pn events. X- and y- axes are set to 
the same scales. We give the fractional r.m.s. variability of each 
lightcurve. 
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Figure 3. Power density spectra constructed from 0.3-10 keV 
pn events for ULXl (black) and ULX2 (grey); the expected noise 
is subtracted. The PDS of ULXl features strong variability at 
low frequencies; we provide the best fit power law model with 7 
= 1.0±0.5 (x^/dof = 3/6). 



(FLAG==0) filter is applied; lastly, ULXS is more heavily 
absorbed than the other two. ULXl is particularly variable, 
ranging in intensity over ~0.0-0.4 count s~^; analysis of the 
ULXl lightcurves in the 0.3-2.5 keV and 2.5-10 keV bands 
showed no evidence for energy dependence in this variability. 

For each source, we constructed 0.3-10 keV PDS that 
were averaged over several intervals consisting of 64 bins 
with a resolution of 100 s; intervals of background flaring 
were excluded. Figure [3] shows the resulting PDS for ULXl 
(black) and ULX2 (grey) , which have comparable inten- 
sities. The PDS are normalized to give the (r.m.s./mean)'^ 
power and the expected noise is subtracted; the PDS are 
binned by to give 56 frequencies per bin. The PDS of ULXl 
(black) is well described by a power law; if the power at fre- 
quency 1/ is P(v), then Piy) tx v~'' . We show the best fit 
power law to the ULXl PDS in Fig. O 7 = 1.0±0.5, with 
X'^/ = 3 for 6 degrees of freedom (dof). Fitting the ULXl 
PDS with zero power yields x^/dof = 26/8, an unaccept- 
able fit; hence the observed power in ULXl is significant. 
The PDS for ULX2 shows no strong evidence for variability, 
even though ULX2 is slightly brighter than ULXl; fitting 
zero power to the ULX2 PDS yields xVdof = 5/6. Hence 
the observed power from ULXl is intrinsic, rather than an 
artefact of the observation. Indeed, the PDS of the combined 
pn+MOSH-MOS2 lightcurves of ULX2 and ULX3 were also 
featureless. These results show that 10,000 source counts are 



sufficient for detecting 7^1 PDS in sources where the r.m.s. 
variability is ~30%, but not for sources with 10% variability. 
The frequency range of our PDS is ~2x lO-'-SxlO"^ 



Hz.lch urazov et al.l (|200lh showed that Cygnus X-1 exhibits 
similar power law noise with 7 ~1 in this frequency range 
during its low state and high state; however, the r.m.s.^ 
power of Cygnus X-1 at ~10~'* Hz is ^^20 times higher in 
the high state than in the low state. This is unsurprising, be- 
cause the high state PDS is characterisd by a 7 ~1 power law 
over the 0.0001-10 Hz range, while its low state is similar, 
except for the interval ~2x 10~^-2x 10"^ Hz where 7 ~0. 
The variability of ULXl is consistent with that of Cygnus 
X-1 in the high state, but not its low state. 

We compared th e 10~^-10 Hz PDS of C ygnus X-1 in the 
high and low states l|Churazov et al.ll200ll ). with the 10"'^- 
10 Hz PDS of other black ho le X-ra y binaries in the high, low 
and steep power law states (|McClin tock & Remillard 20061). 
The high state PDS of Cygnus X-1 has rather more power 
at 10~^ Hz than the high states of other black hole binaries; 
this may be a consequence of its high mass donor. The low 
state PDS of Cygnus X-1 shows variability at 10~^ Hz that 
is as high as, or higher than the PDS of other black holes in 
the low state. It is difficult to distinguish between the high 
and very high states variability alone; we therefore infer that 
ULXl is in a high accretion rate state that may be related 
to the high or very high/ state. 
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Figure 4. Unfolded 0.3—10 keV pn and combined MOS spectra 
of ULXl, modelled witii disc blackbody and power law emission 
components, and suffering photo-electric absorption by material 
in the line of sight. A constant of normalization accounts for the 
differences in pn and MOS responses. 

3.2 Spectral analysis 

The emission spectra of the brightest ULXs are generally 
well described b y disc emission plus a Comptonization com- 
ponent (see e.g. IStobbart et al.ll2006l) . Figure |4] shows the 
unfolded 0.3-10 keV pn and spectrum of ULXl, simulta- 
neously modeled with a disc blackbody and a power law 
component to represent Comptonization; line-of-sight ab- 
sorption is included, and a constant of normalization ac- 
counts for differences in the pn and MOS responses. The 
Comptonized component contributes most of the 0.3-10 keV 
emission; however, the thermal emission dominates the ~l-5 
keV range. The power law component dominates the spec- 
trum below ~ 1 keV; this "soft excess" is sometimes deemed 
unphysical, as the seed photons for the Com ptonization are 
assumed to come from the inn er disc (see e.g. iRoberts et al.l 
l2005l : lGoncai ves fc Soriall2006l ). In this section, we model the 
spectra of ULX1-ULX3 with several alternative models. 



3.2.1 Disc blackbody models 

The emission spectra of distant ULXs are often successfully 
described by disc blackbody models (see e.g. Ilsobe et al.l 
I2OO9I ). However, for these models to be meaningful, they 
must have inner disc temperatures and luminosities that are 
consistent with the black hole mass, A^bh; the temperature 
is governed by M^^^'*, while the lum inosity should not ex- 
ceed ~0.5LEdd (|Mitsuda et al.lll984l ). 

Successful fits to the NGC253 ULX spectra with mean- 
ingful disc blackbody models would indicate that these sys- 
tems were in the canonical high state; disc temperatures 
range over ^0.7- 1.5 keV for most known black h ole binaries 
in the high state (|McClintock fc RemillMdIbOOd '). Therefore 
we simultaneously fitted absorbed disc blackbody models to 
the pn and combined MOS spectra for each ULX; a constant 
of normalisation is included to account for differences in the 
pn and MOS responses. We reject the absorbed disc black- 
body model for ULXl and ULX3. The disc blackbody fit for 
ULX2 is acceptable: N-a = 1.56±0.08x 10^^ atom cm"^ kTin 



= 1.68±0.04 keV, and xV^of = 397/376 (good fit probabil- 
ity, g.f.p. = 0.22). However, the temperature for ULX2 is 
rather higher than known black hole binaries, and is incon- 
sistent with the modelled 0.3-10 keV luminosity of 3x10"^^ 
erg s~^. F-testing shows the blackbody + power law model 
to be a better fit; the probability that this improvement is 
due to chance is 1.4x10"^. 

3.2.2 Slim disc models 

Sli m disc models have a ls o been appl i ed to some ULXs (see 
e.g. lOkai ima et al ]|2006l ) . iKawafuchil (|2003l ) provides a suite 
of tabular models for XSPECjj, including a standard disc, 
thermal and Comptonized slim discs, with and without rel- 
ativistic effects. Each model is described by four parame- 
ters; the mass of the accretor, M, is given in solar masses; 
the accretion rate, M, is normalised to LEdd/ c^, so that 
Eddington accretion corresponds to M ~ 10 (|Kawaguchil 
I2OO3I I: the viscosity parameter, q, is limited to the range 
0.01-1; the normalisation is defined as (10 kpc/d)'^, where d 
is the distance. The distance was fixed to 4 Mpc, following 
iGrimm et all (|2003l ). 

We applied these slim disc models to our ULX targets. 
For ULXl, no thermal model resulted in good fits, meaning 
that Comptonization was required. The best results were 
gleaned from the Comptonized slim disc model with rela- 
tivistic correction: A''h = 1.53±0.06x 10^^ atom cm"^, M 
= 17.5±1.5 Mq, M = 12.8±0.8 LEdd/c^ a = 0.11±0.03, 
and xVdof = 351/324 (g.f.p. = 0.14). For ULX2, all slim 
disc models failed, with g.f.p. <2xl0~''. ULX3 was best 
described by a modified blackbody with relativistic correc- 
tion: A''h = 3.7±0.2xl0^^ atom cm"^, M = 8±7 M©, M = 
10.7±0.7 LEdd/c^ a = 0.36±0.10, xVdof = 94/77 (g.f.p. = 
0.09). F-testing showed that our two-component model has 
a 99.6% probability of providing a significant improvement. 

3.2.3 Self- consistent Comptonization of a disc blackbody 

ISteiner et al.l (|2008l ) have produced a convolution model for 
XSPEC vl2, SIMPL, which describes the Comptonization of 
any input spectrum; the modeling can include up-scattering 
and down-scattering, or upscattering only. It uses only two 
free parameters: the photon index of the resulting power law, 
and the fraction of scattered photons. Fitting our ULXs with 
SIMPL®DISKBB XSPEC emission model will help determine 
whether the observed soft excess is a feature of the ULX 
spectrum or the model. 

Best fit SIMPL models are provided in Table [2] these in- 
clude up-scattering and down-scattering. We reject all SIMPL 
fits to ULXS. We found acceptable fits for ULXl and ULX2; 
however, the emission parameters could not be constrained. 
The best fit to ULXl is presented in Fig.O F-testing tells us 
that we must prefer the two-component models, since they 
provided a lower for the same degrees of freedom in each 
case, and this is always a significant improvement. We there- 
fore infer that the soft excess is a genuine feature of ULX 
spectra; such a feature could be produced by an extended 
corona that has access to low energy photons from the outer 
disc. 

^ http: / /heasarc. gsfc.nasa.gov/docs / xanadu/xspec / models / slimdisk.html 
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Figure 5. Unfolded 0.3—10 keV pn and combined MOS spectra of 
ULXl with the best-fit SIMPL model. A constant of normalization 
accounts for the differences in pn and MOS responses. 



One might expect the SIMPL and two component models 
to give equally acceptable results, as they are both describ- 
ing the Comptonization of a disc blackbody, with no spec- 
ified corona geometry. However, the SIMPL model does not 
include photons outside the energy range of the spectrum; 
hence, the softer photons in the cooler regions of the disc 
are excluded. We further note that the SIMPL model assumes 
that all photons have equal probability of being scattered, 
and that photons of all ene rgies are scattered by the same 
amount (|Steiner et al.ll2008l '). The latter may approximate a 
compact corona, but is unlikely to be true for an extended 
corona, as the photon energies will vary by several orders of 
magnitude. 

We also fitted SIMPL models to our XMM-Newton spec- 
tra for the known black h ole -|- Wolf-Rayet binary ICIO X-1 
fsee lBarnard et al]|2008al . and references within). The best 
fit model yielded xVdof = 763/650 (g.f.p. = 0.0014) and is 
rejected; the two component model yields x'^/dof = 699/650 
(g.f.p — 0.08). Hence we infer the soft excess to be a genuine 
feature of the ICIO X-1 emission spectrum also. 



4 DISCUSSION AND CONCLUSIONS 

ULXl is a factor of ~3 times more variable than ULX2, 
despite the fact that ULX2 is brighter. Furthermore, the 
~0.0002-0.005 Hz PDS of ULXl is weU described by P{u) 
oc ;y-io±0-5; the power at 0.0002 Hz is a factor ~10 
higher than expecte d for the canonical low state for known 
black hole binaries (|McClintock fc RemillardI |2006| ). How- 
ever, the 0.0002 Hz powe r of ULXl is consiste nt with the 
high state of Cygnus X-1 (|Churazov et al.ll200ll '). It is diffi- 
cult to differentiate between the high state and very high 
state (a.k.a. steep power law state) using only the PDS 
iMcClintock fc RemiUardI bOOd ) ; however, the emission of 
ULXl is more like the very high / steep power law state 
than the high state. 

The best studied ULX spectra exhibit are well described 
by thermal -I- Comptonization emission models. They often 
exhibit soft excesses which have been labeled unphysical, as- 
suming that the seed photons are provided by the inner disc 
l|Roberts et al.ll2005l : [Goncalves fc Soriall20odl . e.g.). We have 



tested this assumption using simpl models. We have rejected 
these models for ULXS and ICIO X-1, and find that our two- 
component models are preferred for ULXl and ULX2. We 
therefore infer that these soft excesses are genuine features 
of ULX spectra. 

The disc blackbody contributes 1.3±0.2xl0^^ erg s"^ 
to the unabsorbed 0.3-10 keV flux from our best two- 
component fit to ULXl; a distance of 4 Mpc is assumed, 
but the distance is not well known. This is plausible for a 
black hole with mass ^^20 Mq ; the most massive known stel- 
lar mass black hole, in ICIO X -1, has a mass of ~20-35 Mq 
^Silverman fc Filippenkoll2008l ). 

We propose an explanation for t he soft excess that cru - 
cially involves an extended corona. In lBarnard et all (|2008ar i 
we proposed that the canonical high soft state observed in 
the transient black hole binary systems is caused by the ejec- 
tion of the disc corona as the X-ray luminosity increases by 
several orders of magnitude during outburst. We also pro- 
posed that a steady high accretion rate could maintain a 
stable corona, resulting in a non-thermal high state. This 
theory is supported by observations of the high mass black 
hole binary LMC X- 1, which appears to be in a stable, non- 
thermal high state l| Barnard et al.l l2008al . and references 
within) . 

The "non-standard" fits to the ULXs resemble the emis- 
sion generally seen i n the Galact ic X-ray binaries with neu- 
tron star primaries; IWhite et al.l ([l988 ) found neutron star 
LMXB spectra to be dominated by a Comptonized emis- 
sion component, with an ad ditional blackbody component 

at high luminosities, while IChurch fc Baluciriska-ChurchI 

(|200ir ) found that the blackbody component was present at 
low luminosities also. The thermal and non-thermal emis- 
sion components are spatially distinct in the neutron star 
LMXBs, as demonstrated by the high inclination LMXBs 
known as the dipping sources where the X-ray source is 
experiences increased absorption by material in the outer 
accretion disc on the orbital period; the thermal com- 
ponent is completely and abruptly removed by dipping, 
while the Comptonized component experiences more grad- 



ual changes in absorption (see e.g . White fc Swa nk 1982; 
IChurch fc Balucihska-Churchlll995[|Barnard et alllBoi] ). 

Corona sizes can be estimated for high inclination 
X-ray binaries, from the ingress times of the inten- 
sity dips caused by photo-electric absorption of X-rays 
by cold mat erial in the outer accretion disc (see e.g. 
IChurchI [2OO1I , for details). Measured corona diameters for 
the dipping sources range from 20,000 k m to 700,000 
km with radius increasing with lu minosity l|Churchl I2OOII : 
IChurch fc Balucins ka-Churchll2004l ): the measured coronae 
extend over 10-50% of the accretion disc radius, with this 
fraction increasing with luminosity. Such extended coronae 
would have access to an enormous reservoir of photons -Cl 
keV that could account for the soft excess. 

We note that further, independe nt evidence for ex - 
tended coronae has been discovered bv lSchulz etahl (|2009l ). 
They found broadened emission lines in Chandra observa- 
tions of Cygn us X-2, requiri ng Doppler velocities of '^llOO- 
2700 km s"^ ISc hulz et all 12009} infer from these results 
that the corona of Cygnus X-2 is hot, dense and extended 
up to ~10^° cm. 

If the coronae in ULXs are similarly extended, as our 
results suggest, then the Eddington limit may be relaxed 
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Table 2. Best fit self-consistent Comptonization models to the ULXs. The spectral model used was CONST*wabs*SIMPl{diskbb) with 
CONST accounting for the differences in pn and MOS response. The models include up-scattering and down-scattering. We show the 
absorption (A^h/IO^'' atom cm~^), photon index (F), scattering fraction (/scat), and inner disc temperature (kTin). We then give the 
X'^/dof and good fit probability. The emission parameters are unconstrained for all models. 



Model Afn/lO^O atom cm-2 V /^cat kTi^ xV^of [g.f.p] 



ULXl 6 l.f 0.3 f.O 359/323 [0.08] 

ULX2 16 3.2 0.6 1.3 386/374 [0.32] 

ULXS 24 1.1 0.4 0.8 110/76 [7E-3] 



and many observed ULXs may be explained by stellar mass 
black holes accreting at a consistently high rate, but radi- 
ating at a locally sub-Eddington level. We propose that the 
non-thermal component of the observed ULX spectrum is 
provided by an extended corona that Comptonizes photons 
from the inner regions of the disc: the hot photons from the 
inner disc and also cool photons from further out. The ob- 
served thermal component would then be the portion of the 
disc emission that is unscattered. 

This state may well be related to the canonical very 
high / steep power law state. Indeed, if the high soft state is 
caused by ejection of the corona during transient outburst, 
then the very high state could indicate that the corona was 
regenerating; an extended corona could make the very high 
state appear brighter in X-rays than the high soft state if 
sufficient cool photons were inverse-Compton scattered into 
the observed energy range. It is unfortunate that the highest 
quality X-ray observations of Galactic black hole binaries, 
from RXTE, are insensitive to energies below 2 keV; sensi- 
tive spectra could reveal a soft excess, and extended corona, 
in the very high state also. 

The extended corona scenario means that many ULXs 
could contain stellar mass black holes that exhibit a super- 
Eddington luminosity while being locally sub-Eddington. 
Furthermore, it may be directly tested, as ICIO X-1 is a 
high inclination system that eclipses on a ~35 hour period 
(jPres twich ct al. 2007; Silverman & Filippcnko 2008); if our 
scenario is correct, then we would observe photo-electric ab- 
sorption dips on the orbital period, and could determine the 
size of its corona from the ingress time. 
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